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ABSTRACT 

We present a search for RR Lyrae variable stars from archival observations 
of the Southern Edgeworth-Kuiper Belt Object survey. The survey covers 1675 
square degrees along the ecliptic to a mean depth of 1^=19.5, i.e. a heliocentric 
distance of ~ 50kpc for RR Lyrae stars. The survey reveals 2016 RR Lyrae 
candidates. Follow-up photometric monitoring of a subset of these candidates 
shows ~24% contamination by non-RR Lyrae variables. We derive a map of over- 
density of RR Lyraes in the halo that reveals a series of structures coincident with 
the leading and trailing arms of debris from the Sagittarius dwarf galaxy. One of 
the regions of over-density is found on the trailing arm, 200° from the main body 
of the Sagittarius dwarf at a distance of ~ 45kpc. This distant detection of the 
stellar population of the outer trailing arm of Sagittarius offers a tight constraint 
on the motion of the dwarf galaxy. A distinctly separate region of over-density 
is seen towards the Virgo Over Density. 

Subject headings: Galaxy: halo — Galaxy: structure — stars: variables: other 



1. Introduction 

The stellar halo of our galaxy offers a window to the past. Stars in the halo have 
long orbital periods and inhabit regions of space where the galactic potential is relatively 
smooth and slowly evolving. Consequently, structures in th e halo suffer reduced spat ial and 



kinematic dissipation and may persist for several Gigayears ( iBuUock &: Johnsonll2005l ). Such 
halo structures offer an insight into the process of galaxy formation, arguably one of the most 
important processes yet to be understood by modern astrophysics. 
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One schematic for the process of halo formation is that due to lEggen et aLl (119621 ). who 
found that the stars on the most eccentric orbits were also the most metal-poor. This was 
interpreted as the signature of monolithic collapse of the proto-galactic cloud; the most metal- 
poor, and hence oldest, star s forming during the collapse. This solitary process of galaxy 
formation was challenged by ISearle fc ZinnI (119781 ). who proposed the large dispersion seen 
in the metallicity of globular clusters and their lack of radial abundance gradient indicated 
that the Galaxy was the composite of a large number of sub-galaxy sized pieces. 



Recent observational studies have converged on a combination of scenarios (IFreeman fc Bland-Hawthor 



2002 



f or a review). This is p erha ps most clear l y dem onstrated by the works of iKinman et al. 



(12007). ICaroUo et a 



(|2007l) and iMiceli et al.l (120071 ) who find a distinct outer halo popula- 



tion. iCaroUo et all J2007h shows that whereas the inner halo (galactocentric radius less than 
15kpc) is dominated by highly eccentric, prograde orbits and a metallicity of [Fe/H]~ —1.6, 
the outer halo exterior to this possesses a more u niform distribution o f eccentricities, in- 
cludes highly retrograde orbits and [Fe/H]~ —2.2. iKinman et al.l (120071 ) study of the local 
halo similarly identifies a component that possesses retrograde rotation and streaming mo- 
tions and another that exhibits negligible rotation and less signs of streaming. Kinman et 
al. furthermore show that the horizontal branch morphologies of the two components of the 
halo are consi stent with tho s e seen in the young and old globular clusters, reinforcing the 
conclusions of iLee fc Carneyl (119991 ) that the inner halo is consistent with rapid collapse and 
the outer halo with the ongoing accretion of dwarf galaxies. 

A halo that formed from in situ star formation would possess, at the current epoch 
many dynamical times later, very little substructure. Simulations of galaxy formation in the 
context of ACDM predict that the majority of t he stars in the outer halo formed in progenito r 
dwarf galaxies and were subsequently accreted (lAbadi et al.ll2006l : iBullock fc Johnsonll2005l ). 
with the general conclusion that significant amounts of substructure should survive to this 
time. Howev er, we must point out that ACDM predictions on galaxy scales are contentious 
at this stage. iBell et al.l (120071 ) compare the level of structure evident in the SDSS data with 
simulations and finds that the Milky Way's halo is well matched to a halo built exclusively 
from disrupted satellite galaxies. 

Numerous studies of the spatial distribution of stars in the halo have shown the im- 



The ensuing trail of debris has been found to wrap around the sky ( 



Ibata et al. 


1994. 


Sgr) 


(Belokurov et al. 


2006 



Maiewski et al.ll2003l : llbata et al.ll2003l : lYannv et al.ll2003l : iNewberg et al.ll2002f ). The debris 
of Sgr is arguably the largest contributor yet found to the halo substructure at radii of over 
50kpc. The orbit of Sgr, as traced by its extensive debris tails, has been used as a probe 
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of the mass and shape of the Galaxy 


(Ibata et al. 


2001; 


Helmi & White 


2001 


Johnston et al. 


2005; 


Law et al. 


2005; 


Fellhauer et al. 


2006). 



Another halo structure, the Vir^o Stellar S t ream (VSS) was found as an over-density 



of RR Lyrae stars (RRLs) (j Vivas fc Zinru |2006| ; IZinn et al. 



2004r) at galactocentric radius 



~ 20kpc. Later radial velocity m easurements (jPuffau et al.l 120061) re vealed the grouping 
shared a comrnon sp ace velocity. Newberg et al. ~( |2002 ). Jurid ( 2005 ) and more recently 
Newberg et al.l (120071 ) . find a coincident diffuse over-density of F-type main sequence stars 
spanning some 1000 square degrees of sky called the Virgo Overdensity (VOD). lNewberg et al, 
( 120071 ) proposes that the VSS and VOD are the same structure. On the basis of radial ve- 
locity measurements Newberg et al. concludes it is not associated with a co nfluence of the 
leading and trailing arms of Sgr as proposed by lMartfnez-Delgado et al.l (120071 ) but is possibly 
remnant material from another merger event. 



The Monoceros Stream (INewberg et al.l 120021) is beli eved to be a low surface brightness 
ring encirchng t he Ga laxy at 15-20kpc (llbata et al.l 120031 ). The struc ture was first traced by 
Newberg et al.l (120021 ) using halo turnoff stars and subsequently by iMajewski et al.l (120031 ) 
using M giants. The stream is not thought to be related to Sgr tidal debris or the VSS /VOD, 



but r a ther may result from the disr up tion of the Canis M ajor dwarf (iMartmez-Delgado et al. 



20071 ; ICasetti-Dinescu et al.l l2006l ). I Conn et al.l (120071 ) reports that interpretation of the 



Monoceros stream as a warp or flare in the galactic disk does not account for the observed 
stellar densities. P ossibl y related to the Monoceros Stream are the Tri-And over-density 

20kpc and a yet more distant stellar structure at 



Rocha-Pinto et al. 


2004 


) at 


Martin et al. 


2007 


). 



28kpc 



Most recently, iBelokurov et al.l (120071 ) reports the detection of a significant stellar over- 
density that covers some 80° of galactic latitude, dubbed the Hercules-Aquila Cloud. This 
over-density lies at a distance of between 10-20kpc and straddles the Galactic disk offset by 
6-15kpc from the Galactic centre. Spectroscopy of probable cloud members shows that the 
structure possesses a radial velocity offset of ~180km/s from the halo and thick disk. 

The features discussed above are a small fraction of the number that could be ex- 



pected from simulations of ACDM hierarchical structure formation (iBuUock et al.l 12001 



Bullock fc Johnsonll2005l ). The search for relic substructure could help illu minate the out- 

standing 'missing satellite problem' in galaxy formation in ACMD cosmology (IFreeman fc Bland-Hawthorn 
2OO2I ). Here, models predict of order 500 small dark matter halos should surround the Milky 
Way today, in contrast to the ~ 40 known dwarf galaxies in the Local Group (IWilkinson 
20071 ). If the level of substructure does not match that predicted from a large population of 
dwarf galaxies this may prompt revisions to model assumptions. 
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2. RR Lyrae Variables as Probes of the Halo 



RRL variables have a long history as probes of the Galactic halo (see lWetterer fc McGraw 

(Il996l ) for a historical review). They are stars on the core He-burning horizontal branch that 



fall within the temperature range of the instability strip. They are easily identified on the 
basis of their color and variability. In addition, their absolute magnitudes are known to 
better than 0.1 magnitudes, making them standard candles for the halo. They represent 
an old stellar population and so by virtue of their age can probe disruption events that 
occurred at early times. The abundance of RRLs is dependent on the morphology of the 
horizontal branch of the progenitor population. Local dwarf spheroidals exhibit abundant 
RRL populations, hence these variable stars should be good tracers of accretion events. 

In recent years several sky surveys have discovered large numbers of RR Lyrae vari- 
ables. The survey area, li mits and completenes s of these surveys are summarised in Table [TJ 
The QUEST RRL survey (I Vivas &: ZinnI 120061 ) proposes six significant substructures; three 
associated with the Monoceros Stream, and others associated with Pal5, VSS, and the Sagit- 



tarius Stre am. Searche s for R RLs in the SDSS are presente d bv ISesad (20071). Ilvezic et al. 



(120041 ) and Ilvezic et al.l (120001 ) from multiple imaging and by Ilvezic et al.l (120051 ) from colors 
alone. ISesarl (120071 ) reviews the SDSS work to date, proposing 13 new over densities. 



The QUEST a nd SDSS RRL surveys cover equatorial regions. The LONEOS-1 Survey 
( iMiceli et al.l l2007l ) traces the ecliptic but to a moderate depth ( reaching distances of up 
to 30kpc). With between 28-50 epochs per star lMiceli et al.l (120071 ) are able to identify the 
Oosterhoff type of each RRL. They conclude that the radial di stribution of O o sterho ff type 
I and II (Ool and OoII) RRLs are systematically different. iLee fc Carneyl (119991 ) found 
globular clusters that are dominated by Ool RRLs show li tt le ne t angular momentum as 
expected from a formation scenario such as ISearle fc ZinnI (Il978l) whereas Pol l globular 
clusters show net p rograde rotation as would result from the lEggen et al. ( 1962 ) scenario. 
Miceli et al.l ( 120071 ) propose that the presence of distinct Ool and OoII populations reveals 
that both the monolithic collapse and accretion scenarios are important for constructing 
the halo. This con clusion parallels the two kinematic components found in the halo by 



CaroUo et al 



m 



Miceli et al 



{ 20071). The presence of substructure in their sample of RRLs is not addressed 



(!2007f ). 



3. The Southern Edgeworth-Kuiper Belt Object Survey 



The present study is based upon t he archived observa tions of the Southern Edgeworth- 
Kuiper Belt Object (SEKBO) survey (IMoody et al.l 120031 ). SEKBO was a multi-epoch, two 
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color imaging survey of a 10 degree band centred on the ecliptic. Observations were obtained 
over the three years from December 1999 on the Mount Stromlo Observatory 50-inch Great 
Melbourne Telescope that was previously used for the MACHO microlensing survey. The 
telescope featured a dual-band imaging system in which a dichroic allowed for simultaneous 
exposures in two broad bands on separate 4kx4k pixel CCD mosaics. The MACHO blue 
filter, hereafter Bm, spans 455-590nm (wavelengths at 50% normalised response) and the 



red, Rm, from 615-775nm (iBessell fc Germany! 1 19991 ). These wide, non-standard filters most 



closely match Johnson V and R. 

SEKBO was designed to search for Trans-Neptunian Objects from a set of three expo- 
sures of a given field, the first and second epochs separated by ~4 hours and a subsequent 
third epoch 1-7 days later. As the telescope operated in an automated manner with limited 
environmental monitoring some observations were obtained under conditions of insufficient 
quality. These would break an observing sequence for a field and require reacquisition of 
the field at a later date. Figure [1] shows the distribution of time intervals between repeat 
observations of the fields in our sample. The exposure time was 300 seconds. Under rea- 
sonable photometric conditions, this results in a useful magnitude range for photometry of 
14 < < 20. Figure [2] shows the distribution of seeing derived from the SEKBO images. 
The seeing is typically 2" with a tail out to > 5" . 

There are 5212 field centres with two or more epochs that form the basis of our sample. 
The histogram of the number of epochs per field in the sample is shown in Figure [31 Most 
fields have 2 or 3 observations but there is a tail out to 10 observations. Figure S] shows the 
sky coverage of the survey. Note that we have excluded regions with E{B — V) > 0.10 as 
well as galactic latitudes b < 10° to avoid fields with excessive crowding. The total effective 
area of the survey (having corrected for inter-CCD gaps and bad columns, see ^6]) is 1675 
square degrees. 



4. Photometric Reduction 

In this section we outline the steps taken to reduce the raw frames to calibrated pho- 
tometry. Each of the 17578 exposures considered consists of 16 images (2 colors x 4 CCDs 
X 2 amplifiers) which are stored on the Australian National University's Supercomputing 
facility mass data store. The data set is 2.3Tb in size. 

The data in the archive has been processed to remove bias and fiatfield structure using 
the data reduction pipeline of the MACHO system. This system maintained a moving 
average bias and fiatfield that was regularly checked for changes in the optical/detector 
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system dAlcock et aDll999h . 



We first determine a world co-ordinate system (WCS) for eacli image based upon tlie 
nominal pointing of t fie telescope as recorded in the image FITS header. UCAC2 stars 



(jZacharias et al.ll2004l ) in common to the image are used to derive the WCS. The UCAC2 
provides between 40 and several hundred stars per CCD. Our derived positional uncertain- 
ties are better than 50 milliarcsecs for bright stars. A ZP]n|i] coordinate representation is 
used. This representation utilizes a radial term to characterise distortions across the image 
and a simple linear transformation to remove the effects of rotation, scale, and any shear 
caused by refraction. Linearising the telescope system in this way provides robust coordinate 
transformations across multiple images of a field with no difficulties in CCD corners or in 
images with few astrometric standards. 

Differences in the fabrication and electronics associated with each CCD/amplifier pair 
result in a non-uniform gain across the red and blue mosaic. Gain matching between each 
CCDs two amplifiers was achieved by matching the median sky background level. The four 
CCDs were then mosaiced together using the SWarp packagcl. The median sky brightness 
was then used to scale the four CCDs of each passband to the same gain level. 



4.1. Photometry and Object Classification 

We base our object catalogue upon a template image that is created for each field by 
co-adding the individual Ru observations. In the case of three or more images we use SWarp 
to form a median image to robustly reject bad pixels and cosmic-rays. In the case of two 
images for a field we apply an algorithm that examines corresponding pairs of pixels in the 
registered image stack and removes high values from pairs with large variance. The footprint 
of the template image is only that area on the sky that is represented in all images of the 
field. 

The template image is used solely for object detection. The combination of images taken 
under different seeing would render the extraction of photometry from the template prob- 
lematic due to the distorted and non-photometric PSF shape that results from combination. 
The template image does have several advantages for object detection, however, including 
greater depth and reduced numbers of false objects due to the exclusion of cosmic-rays 



^Zenithal/azimuthal polynomial projection see ICalabretta et al.l (|2004J l 



^http : //terapix . iap . fr/. lBertinI (|2005l ) 
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and bad pixels. Object detection within the template image was made using the cdoPhotlfl 
package. 

Photometry was then performed on the individual images that comprised the template. 
This photometry was performed in three modes: in fixed position aperture mode where the 
template coordinates are used for the centre of the aperture; centroid aperture mode where 
the barycentre of the object is used for the centre of the aperture; and PSF-fitting mode 
where a PSF is constructed for each object. Aperture corrections were applie d to corr e ct the 



magnitude defined through a 3 pixel radius aperture to the aperture used by iLandoltl (1l992l ) 
for standard star photometry (discussed below). 

The PSF-fitting mode is the basis of our object classification. The gaussian widths 
in the x and y directions give a clear discriminant as to the nature of the object. Stellar 
objects lie clustered around a particular x-width, y-width. Galaxies exhibit more extended 
dimensions and cosmic-rays a smaller extent. To make a quantitative judgement as to the 
stellar nature of the objects in our data we construct a stellarity index that expresses the 
distance an object lies in gaussian x-width, y-width space from the stellar locus. 



4.2. Standardised Instrumental Photometry 

Raw magnitudes from the gain matched frames were then placed on a standard instru- 
mental system. Due to non-photometric conditions and atmospheric extinction there were 
zero-point offsets between magnitudes in each observation and that obtained from the detec- 
tor system at the zenit h under photometric c onditions (hereafter the standard instrumental 



system). The 2MASS (ISkrutskie et al.ll2006l ) J and K magnitudes allow this offset to be 



calculated. From images obtained in photometric conditions we first derived a linear fit: 

RM,i -J = -2.5 log /^j^, + 2.51ogtexp -J = aiJ-K) + ZPr^. (1) 

where f^^ . is the flux in MACHO R for the ith observation, texp = 300 seconds and ZP^^ . 
is a zeropoint defined below. We determine a = 1.100 over a narrow color range (0.3 < 
J-K < 0.6). 

For each individual image in the data set, a fit was then performed on the order of 100 
stars in each field with signal-to-noise ratio > 20 (J < 15 and K < 14.5) within the color 
range above. We define our standard instrumental system such that a Rm = 22 star provides 



■^http : //www . vuw . ac . nz/ staff /michaeljreid/ software . html 
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1 ADU s-\ Hence: 

Rm,! — ZPr^. = 22 (2) 

where ZP^^. is the zeropoint required to draw the observation onto the standard instru- 
mental system (ZPr^ . is the extrapolation of the linear relation of Equation [T] to J—K=0). 

The effective wavelength of the Bu filter is sufficiently blue of the 2MASS bands that 
the above procedure can not be used as it would introduce a significant interstellar ex tinction 



deperi dance. To determine the Bm zeropoint we used the V—R color of sources from [Stetson 



(I2OOOI ) and their corresponding Bm—Rm colors. Note that because the blue and red exposures 
are obtained simultaneously there is a single offset between Bu and Ru (with a minor airmass 
dependence, see below) and this is found by the requirement that 

V-R = = iBM- ZPbJ - {Rm - ZPrJ (3) 



Since the MACHO filters are well matched to the standard V and R only a linear 
regression is required. However, due to the chromatic effect of atmospheric extinction, the 
difference in zeropoints will be a function of airmass, X. From the observed Landolt fields 
we derive, 

ZPbm - ZPr^, = (0.316 ± 0.06) - (0.08 ± 0.02) (X - 1) (4) 
and scale the flux of the blue images to provide 1 ADU s^^ for a Byi = 22 star. 

We reject those observations displaying zeropoint magnitudes that show more than 0.75 
magnitudes of extinction or sky brightness that are more than 1 magnitude brighter than 
the median. We also remove those observations with seeing in excess of 3.6" to reduce 
the problems introduced into the RRL detection technique by objects that are blended in 
poorer seeing frames and resolved in nominal seeing frames. This leaves 3692 field centres 
with acceptable data. 



4.3. Transformation to V and R 



A number of fields in the dataset contain standard stars observed by IStetsonI (I2OOOI ). 
This allows us to define the transformation equ ations to take our standard instrumental 
photometric system to that of Landolt V and R (lLandoltlll992l ). Given the close matching 
of the MACHO bandpasses to the V and R filters only a linear transformation in color is 
required. The transformation is of the form: 



V = BM + aiBM-RM)+P 



(5) 
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for constants a and /5, and where -Bm and -Rm have been both zero-point offset and aperture 
corrected. Observations of a large sample of standard stars allow us to define the following 
conversions to Landolt V and R magnitudes 



and. 



V = B^- 0.12 (5m - i?M) + 0.19 



R = Rm + 0.15 {Bm - Rm) + 0.34 



(6) 



(7) 



Figure [5] s hows the differe nce between the magnitudes in V and R derived from Equation 
inland those of IStetsoru (120001 ) from the 160 stars in common. This figure shows an offset of 

= 0.01 ± 0.03 and Ai? = 0.01 ± 0.03. It should be noted that the photometry that 
contributes to Figure [5] was obtained under hetrogeneous seeing, atmospheric extinction and 
photometric conditions over the course of three years giving us confidence in our global 
zeropoints to the level of ±0.03 magnitudes. 



5. RR Lyrae Candidate Selection 

We have implemented a probabalistic approach to the selection of RRL candidates. 
We first remove those objects with shape and color parameters grossly unlike the stellar 
population. Selection is then made based on color and variability. In order to quantify the 
significance of variation we need to first understand the photometric uncertainties. 

Photometric uncertainties were charactersised using the observed standard deviations 
of the stellar objects from the multiple epochs. Figure [6] shows the distribution of standard 
deviation versus magnitude for objects in a particular survey field. For most stars the stan- 
dard deviations are indicative of the photometric uncertainties for the observation. However, 
variables stars or objects with bad photometry on some epochs will lie above the main locus. 

To determine the mean photometric uncertainty as a function of magnitude we fit a 
sky-limited Poisson distribution: 

a,Hot = 10^'-^^"'+"'^-^^ + m floor (8) 

for magnitude m, where rricut is iteratively determined so that the standard deviation of half 
the objects lie above the fitted line and half below. The constant m floor is the component of 
the photometric uncertainty that is not related to the photon statistics. This is derived from 
the brighter stars up to two magnitudes fainter than the saturation limit. The photometric 
uncertainties vary from field to field. The typical uncertainty for the brighter stars is 0.02 
magnitudes rising to 0.2 magnitudes at ~ 19.5. 
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With an understanding of our photometric uncertainties we then construct a RRL index, 

£: 

= ^^Color + ^^RRL — ^NotVariable (9) 

where Ccoior quantifies the probabihty that a star has colors commensurate with an RRL; 
^RRL is the log likelihood that the variation in color exhibited by the object matches that 
expected from RRLs; and C-Notvariabie is the log likelihood that the variation can be explained 
by photometric uncertainty alone. The RRL index, £, is a figure of merit of how much like 
an RRL each object is. It is not based on a deep statistical foundation, but rather the 
pragmatic need to distinguish RRLs from other stars with a minimum of observations and 
with minimal manual intervention. 

We now discuss the above three quantities in more detail. To construct Ccoior we first 
derived the de-reddened color distribution from ~ 19000 LMC RRLs from the MACHO 



project (using a mean reddening of Fj(V—R) =0.08 towards the LMC (lAlcock et al.ll2004l )). 
This distribution can be approximated by a gaussian with {V — R)o = 0.137 and ay-R = 0.08 
mag. Given this distribution the probability density of a star having an RRL-like color is: 



2 ^ D ^ 



Ccoioriy - R\i = j=e -v^n ' (10) 

where the dereddened V — R color of the i-th observati on, {V — -R)n} , is de rived from 



Equation [6] and dereddened according to the dust maps of ISchlegel et al.l (119981 ) (we adopt 
E(V-i?)=0.608E(S--l/) and Ay=3.315E(5-1/)). 

i^RRL compares the observed variation for an object to the variations seen from a sample 
of RRLs. Densely sampled (Nobs ~ 1000) -Bm and Rm light curves were obtained for a 
random sample of MACHO LMC RRLs. For each point in the light curves we calculate 
Ai?M = Bm,! — {Bm) and ARm = Rm,! — (-Rm)- These are plotted in Figure [71 The line 
ABm = 1.32A_Rm shows a classic signature of pulsation, namely, that the amplitude of 
pulsation increases with decreasing wavelength. To compare these points to the observed 
variation in our objects we determine the fraction of MACHO points enclosed within an error 
ellipse centred on {ABm,ARm) in Figure [7] for each epoch. CNotVar measures the probability 
that a non-variable star of given magnitudes -Bm, -Rm and photometric uncertainties as, ctr 
could have the observed variation. CNotVar for the i-th observation is expressed as follows: 

CNotVar,^ = 6 '^R > € > (11) 

iTXaBUR 

The combined RRL index, £, for a typical field is shown in Figure [81 Of the 3928 objects 
identified in this field only 523 objects pass the shape and color restrictions and remain in 
this figure. The bulk of stars without RRL-like variability and color have indices less than 
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zero. Two objects in the field in question show high indices. We identify RRL candidates 
above the dashed hne in Figure [3 This line is derived from the distribution of 20 QUEST 
RRLs in our survey area (all are RRa6s). Of these RRLs 8 reside in the non- variable locus 
and 12 present £ > 40. The line slopes upwards with V to avoid the flaring of the stellar 
locus at fainter magnitudes. Candidate RRLs were then visually inspected to remove objects 
associated with obvious image flaws (ghosting, bad columns, proximity to bright stars etc.). 
We also limit candidate colors to V—R < 0.3 as discussed in §7. This resulted in 2016 RRL 
candidates. 



6. Survey Completeness 

The survey completeness is the fraction of RRLs recovered by our detection technique 
from those available in the fleld. A detailed characterisation of the completeness is necessary 
for the analysis to follow. The completeness depends critically on magnitude since our ability 
to recover RRLs decreases as we delve fainter due to increasing photometric uncertainty. In 
addition, the limiting magnitude varies across our surveyed flelds. Therefore it was necessary 
to calculate the completeness and a function of magnitude for each fleld. 

To measure the completeness we created artiflcial RRL light curves drawn from the 
period, amplitude and color distributions of the MACHO LMC RRL sample discussed in 
Section [51 Since RRLa6 and RRLc stars have d ifferent p eriod and amplitude distributions 



we considered the two populations separately. iLayderu (119981 ) template light curves were 
used to construct the RRLaft light curves and sine curves were used for the RRLc stars. The 
initial phase of the simulated light curve was chosen at random. Subsequent samples of the 
light curve were made at the intervals of the observations for the fleld in question to account 
for the phase coverage of the fleld in question. 

To these idealised light curves we added gaussian-distributed photometric uncertainties 
as described by Equation [SI To measure the completeness as a function of magnitude we 
generated artiflcial light curves in 0.5 magnitude bins from 14< Rm <21. The artiflcial data 
were passed through the selection pipeline and the recovered fraction formed the complete- 
ness estimate for that fleld and magnitude bin. The artiflcial data however, do not suffer 
from the imposition of bad columns, occasional dead ampliflers and inter-CCD gaps. To cor- 
rectly account for these features we determine the effective area of each fleld. The effective 
area is the area of the fleld that has science data from each exposure. On average, this is 
89% of the maximum 0.51deg^. 



Typical completeness proflles are shown in Figure M For RRLa6s our survey averages 
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~ 60% completeness for V < 18.5. This falls to a mean completeness of 25% hj V = 19.5, 
although a handful of fields extend to V = 19.5 due to a combination of low sky brightness 
and exceptionally good seeing. The completeness for RRLcs is systematically lower than that 
of the RRLafe due to their lower amplitudes. Not surprisingly, the completeness increases 
with more observations of the field. Figure [10] shows a map of the completeness over the 
survey fields. Fields with anomalous low completeness typically have bad photometry or are 
derived from observations that are separated by a few minutes and are thus insensitive to 
the timescale of RRL variation. 

The overlap between our survey and QUEST provides a check on our determined com- 
pleteness. As stated above, of the 20 QUEST RRLafo stars in our survey area, 12 are 
recovered. This is a magnitude-averaged completeness of 60%. This compares reassuringly 
well with the average completeness expected from these fields of 55%. In the case of RRLc 
stars the numbers are too small to draw any conclusions; one was recovered out of three 
expected. 

6.1. Magnitude Offsets 

The limited number of observations and the asymmetric shape of the RRLa6 light curve 
introduces a bias in the mean magnitude of selected candidates. This can be understood quite 
simply. Due to the shape of the RRLafe light curve less time is spent near maximum light. 
Yet an RRLafo has highest likelihood of selection (highest C) when observed at maximum 
and then again at minimum. This introduces a net bias towards brighter mean magnitudes 
amongst the RRLs selected as candidates than would be derived from well sampled light 
curves of the same stars. For bright RRLs the offset is minimal - these objects are detected 
with high significance from variations that were only a small fraction of the full amplitude 
due to the low photometric uncertainties. As we proceed to fainter candidates, the fraction of 
the full amplitude that must be exhibited for the object to stand out as a variable increases 
as the photometric uncertainty rises. 

Since the photometric uncertainty as a function of magnitude varies from field to field 
we construct the mean magnitude offset for each individual candidate. A series of artificial 
RRLafo light curves was generated at the magnitude of the candidate as discussed above. 
We then determined the offset between the mean magnitude and the mean recovered from 
the observational epochs. The mean magnitude offset for the ensemble of RRL candidates is 
shown in Figure [TTl RRLc stars do not exhibit any significant bias due to their symmetric 
light curves. Since we do not know apriori whether a candidate is ab or c type we weight 
the mean magnitude offset by the fraction of each Bailey type to the RRL population. We 
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use a ratio of RRa6 to RRc of 91:9 ratio (ISmithlll995l ). although this fraction varies with 
the Oosterhoff type of the system. For example, in the case of the Draco dSph (Ool) RRc 
contribute 11% of the 146 known RR Ls (IBonanos et al.ll2004f) whereas in Bootes dSph (Ooll) 
they comprise 7 of 15 known RRLs (jsiegel 2006 ). An OoII halo with 50% RRc would mean 
that we have systematically underestimated the brightness of each object by at most 0.05 
magnitudes towards our faint limit. 



7. Survey Contamination 

In order to address the level of contamination by non-RRL objects we performed follow- 
up observations of a sample of 66 RRL candidates with the Wide-Field Imager on the RSAA 
40-inch telescope at Siding Spring Observatory. The selected candidates spanned a range of 
magnitudes and V—R color. Light curves of b e tween 8 and 16 points per star were obtained 



and examined for periodicity (see iPrior et al.l (120071 )). The break down of this sample as a 



function of color is shown in Figure [T21 Those that could be phased were either RRLs or 
eclipsing binaries. A small number of objects were clearly variable yet unclassifiable and a 
similar number showed no significant variation over the 4-6 nights of observation. 

In the analysis to follow, we apply a color cut to the RRL candidates of V—R<0.3. This 
cut represents a 2 sigma limit of the MACHO LMC RRL color distribution. With such a 
color cut we find 16 non-RRL amongst 66 candidates (i.e. a contamination rate of 24±12%). 



8. Properties of the RR Lyrae Candidates 



The selection process discussed above results in 2016 RRL candidates. Figure [131 shows 
the color-magnitude diagram for the candidate RRLs relative to the general halo background. 
The halo color-magnitude diagram shows a steep drop off in density for stars of V—R <0.3; 
this is the color of the halo main sequence turn off. The distribution of candidates is bound 
to the red by our color limit at V—R =0.3 and to the blue is consistent with the distribution 
of MACHO LMC RRLs. 

We adopt an absolute magnitude of My = 0.56 for the RRL population. This value 
is based upon the mean absolute magnitude of My=0.59±0.03 for an [Fe/H]=-1.5 provided 



by ICacciari &: Clementinil (120031 ) who reviewed the absolute magnitude determinations for 

liques. The mean metallicity of a sample of halo RRLs from 



Suntzeff et al. 


( 


1991 


) is \ 


of 


Cacciari & Clementini 
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assumed absolute magnitude is equivalent to that used in the study of I Vivas &: ZinnI (120061 ) . 
Vivas & Zinn (ibid) find that the effects of evolution away from the zero-age horizontal 
branch and dispersion in metallicity should lead to a dispersion in My of 0.13 magnitudes, 
or in other words a 7% dispersion in distance. 



8.1. Spatial and Radial Distribution of RR Lyraes 

In Figure [H] we present the radial distribution of RRL candidates. We collapse the 
distribution in ecliptic latitude, the natural system for our sample. A series of over densities 
can be seen, in particular, at ecliptic latitude 120°, 190° and 320°. These features are 
however, overshadowed by the spatial non- uniformity of the survey and the general density 
falloff away from the Galactic centre. In the next two sections we characterise the general 
candidate background in order to subtract it and reveal the significance of these features. 

Figure [TBI shows the position of these over densities on the plane of the sky in the context 
of known significant substructure. We now discuss these features in more detail: 



8.2. Spatial Density of RRL Candidates 

It is clear from the substructure evident in Figure [T3] that the halo does not exhibit 
a smooth radial density profile. In order to characterise the significance of the observed 
substructure we determined the general halo background upon which the substructures are 
overlaid. 

Note that our sample includes significant contamination by non-RRL objects. It is 
unlikely that the contaminants share the spatial distribution or absolute magnitude of the 
RRL population. In the analysis that follows we treat each object as an RRL. The effect this 
has on the derived spatial density of candidates depends on the nature of the contaminants. 
For example, main-sequence eclipsing binaries of significantly lower absolute magnitudes 
than RRLs will appear over a range of magnitudes out to the edge of the disk along a line 
of sight. Assuming RRL absolute magnitudes for these objects will act to enhance the space 
density at large distances. The goal of this exercise is not however to derive the spatial 
density of RRLs but rather to derive a model for the background halo. 

There are two key models discussed in the literature for the radial distribution of matter 
in the halo. General ellipsoidal isopycnic contours for the halo can be expressed as: 

- + ^ + - = 1 (12) 
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The spherical model {a = b = c) is the most basic, in which the density varies as a function 
of galactocentric distance as a power law, 



p{R) = po{R/Ro)\ 
where po is the local space density of RRLs and Rq is the Solar radius. 



(13) 



The case of a spheroidal halo has a = b and flattening described by g = c/a. Numerous 
observational studies have found that the i nner < 20kpc of the halo is significantly fla ttened 



while the outer regions appear spherical Jchiba fcBgCT] I2OO0I: 



rest on et al.l Il99ll ) . The 



c/a 



spatially variable flattening ratio introduced by iPreston et al.l (Il99ll ) expresses isopycnic 
contours as ellipsoids of revolution. Such surfaces vary as a function of the semi-major axis 
a, 

(c/a)o + [1 - (c/a)o](a/a„), a < a„, ^^^^ 
1, a> a„, 

where a„ = 20kpc and canonically, (c/a)o = 0.5. In this case, density varies as a function of 
semi-major axis as: 

p{a) = poia/RoT, (15) 



Numerous studies of the RRL distribution have found that the observed distribution is 
best fit by variable flattening model. However, the difference in t he quality of the res ultant 
fit between the sph erical and variably fiattened models is of little (IVivas fc Zinnll2006l ) or no 
(IMiceli et al.l 120071 ) statistical significance. In the following we only consider the case of a 
spherical halo. 



The calculation of the RRL space density follows the technique developed by lWetterer fc McGraw 



(119961 ). The following equation describes the density as a function of galactocentric distance: 

1 dN 



p{R) 



AnR^fiR) dR' 



(16) 



where f{R) is the fraction of the total halo volume at R that is sampled by the survey 
{f{R) is analogous to a solid angle) and N is the number of RRL as a function of distance. 
Whereas the solid angle is constant as one looks through the halo, f{R) varies as a function 
of Galactocentric radius and hence must be calculated numerically for each field. This is 
because the volume of each field is a pyramid subtending the field's effective area (see §6) 
at infinity and an apex offset from the Galactic centre by Rq. Different lines of sight sample 
different regions of the Galactocentric sphere. This is demonstrated for two lines of sight in 
Figure [T6l One field is of higher galactic latitude and samples more of the inner galactic halo. 
To generate these profiles we calculate R and the volume of each slice of the pyramid for a 
given r + dr. The volumes are summed into O.lkpc bins before being divided by the volume 
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of the sphere at the radius R of the bin. As R tends to large values f{R) approaches the 
asymptotic limit of (effective area)/47r, corresponding to the constant solid angle as would 
occur if Rq = 0. 

To account for the effects of completeness we mulitply each f{R) by the Monte-Carlo 
derived RRLafo completeness profile for the field as a function of galactocentric radius. This 
gives the effective volume of the halo that is sampled by each field (the right-hand plot of 
Figure [TBI) . The total volume of our sample is then simply the sum of f{R) for the 3692 
individual fields. This is shown in Figure [T71 

We then use Equation [16] to calculate the local space density for each candidate. To 
distill the space density of our sample we then bin space densities of the 2016 RRL candidates 
in bins of 0.02 dex in \og{R). The resulting average density distribution is shown in Figure 
[T51 The error bars in Figure [TS] represent the standard deviation of the points in each bin. 
Since 50 to 30 objects reside in a typical bin the large standard deviation is indicative of 
intrinsic azimuthal variability in the space density such as would result from substructure. 

As stated above, our goal in the examination of the average space density is to extract 
the mean background upon which the substructure lies. A linear least squares fit to points 
between 10 < i? < 45kpc yields n = -2.48 ±0.09 and po = (6.22±0.12)kpc-3 (as introduced 
in Equation [T2l) . Beyond R=45kpc there is a more rapid drop in the number of candidate 
RRLs. For the 122 candidates with R > 45kpc we find a = —5.0 ± 0.3 power-law. 
We plan to follow-up R>45kpc candidates to see if the RRLs amongst them match this 
increased radial drop-off. However, as stated before, since the candidates contain significant 
contamination by non-RRLs, we do not interpret these results further at present but rather 
use them as the basis of our detection of spatial departures from the mean background. 

8.3. Halo Substructure 

In this section we quantify the significance of substructure in the halo. We firstly 
construct a series of artificial halos based upon the spatial distribution function discussed 
above. This is complicated in the case of our survey by non-uniform spatial, temporal and 
radial coverage, however all three are well characterised for each field. To form a simulated 
halo we first randomly select a position in the sky from the area covered by our survey 
fields. We then use the spatial density power-law indices discussed above to randomly assign 
a distance and then use the completeness information for the field in which it lies to see 
whether we would recover such a candidate. This is repeated until the total number of 
candidates observed is recovered. An example of a simulated halo is shown in Figure [T31 
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We next perform a Voronoi tessellation (jVoronoil 119071 : iRamella et al.ll200ll ) of the sim- 
ulated halo distribution in (heliocentric distance, ecliptic longitude) space. The Voronoi 
tessellation associates with each point a polygon such that every point inside the polygon 
is closest to the point in question. It enables us to define the local density of points in a 
non-parametric way as simply the reciprocal of the area of the tesserae. We repeat this for 
the entire set of simulated halos and form a map of the mean density of simulated candidates 
in (distance, ecliptic longitude) space. Along with the mean we determine a map of the stan- 
dard deviation in the simulated density. From the Voronoi tessellation of the observed RRL 
candidate distribution we subtract the mean density of simulated candidates and divide by 
the standard deviation to obtain a map of the significance of spatial structure in the halo as 
seen in Figure! 



8.3.1. Sagittarius Dwarf Debris Streams 



In the interpretation of Figure [11] it is useful to consider the position of the orbital 
plane of the Sagittarius dwarf galaxy that features prominently in the outer halo (as seen 
in Figure [T5|) . Our survey cut s obliquely through the plane of the Sagittarius dwarf (Sgr; as 
defined by llbata et al.l (120011 )) at an angle of 16° with an opening angle of 10° (see Figure 
Hn]). We obtain maximum signal from structures in the Sagittarius Stream at angles adjacent 
to the line of nodes. The structures from RA=45° to 130° and 220° to 330° lie in this range. 
In Figure [T9] we show a view of the Sagittarius tidal arms compiled from t he observations 



of 



of iMaiewski et al.l (120031 . M giants) , iBelokurov et al.l (120061 . F dwarfs) and iNewberg et al. 



(I2OO3I. I2OO7I. Blue horizon t al branch s t ars an d F d warfs) together with t he model predictions 



Helmi k Whitd (l200lh . lLaw et all (l2005h . and iFellhauer et al.l (120061 ). 



The models of Law et al 



M giants (IMajewski et al.ll2003l ). iHelmil (120041 ) finds that the data strongly su ggests that 



20051) ari d iHelmil (120041 ) aim to match the observations of 



the da rk matter halo of the Milky Way is prolate ((1.25 < {q = a/c) < 1.5). iLaw et al. 
(I2OO5I ) similarly require the presence of a significantly prolate halo (g = 1.25) but could not 
establish a model that reproduces both the observed orb i tal po le precession and the distance 
and radial velocity of the leading arm. iJohnston et al.l (120051 ) conclude that it is necessary 
to consider an evolution in the orbital parameters of Sgr over the timeframe of the last few 
orbits (~2-3 Gyr) and propose that dynamical friction is the most favourable mechanism to 
bring about this change. 



Fellhauer et al.l (120061 ) models the Sgr system in the light of the IBelokurov et al.l (l2006l ) 



result. In particular, Fellhauer et al. seeks to interpret the observed bifurcation of the Sgr 
stream. This work is based upon the distances derived to the F dwarf population. As pointed 
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out by iNewberg et al.l (120031 ) a 13% discrepancy in distance exists between the distances 
derived from M giants and those from F dwarfs. The M giant distance scale has proven to 
be systematically smaller than the F dwarf scale and it is unclear which is more accurate. 
Fellhauer et al. finds that the observed bifurcation is only possible in a near-spherical halo 
and that it is not possible to match the distances and radial velocities of the leading arm 
in their closest approach to the Solar circle. We must bear in mind these caveats to the 
modelling of the debris from Sgr as we now discuss the over densities found in the present 
study with reference to the orbit of Sgr. 

The Sgr leading arm is potentially seen in the over density located at RA=115° and 
distance of 15kpc (see Figure [19]). This region is a confl uence of the Monoceros Ring and 



Sagit tarius Stream both located at 15-20kpc in this region fllbata et al.ll2003l : iBelokurov et al. 



20061 . respectively). The Sgr leading arm then passes above our survey area until it is again 
intercepted at RA=220°, d=50kpc. This region is the southern limb of the S g;r leading arm 
(see F igure [T5|) the northern edge of this structure is seen in the QUEST survey (jVivas fc Zinn 
20061 . their Figure 12). The main body of Sgr lies in the region of the Galactic plane avoided 
by our survey, however, we see over densities coincident with the trailing arm from RA=287° 
to 330°, d=25kpc. The trailing arm then passes underneath our survey area before potentially 
emerging at RA=130°, d=45kpc. 

The over-density at (RA=130°, d=45kpc) is particularly interesting. It is separated by 
~ 200° from the main body of Sgr alo i ig the trailing arm. Other than the over-density of A- 
type stars detected by lNewberg et al.l (120031 ) . this over-density is the most distant and oldest 
portion of Sgr debris yet found and offers a critical constraint on the orbit of Sgr and hence 
the potent ial of the Galaxy in wh ich it moves. As we ca n see from the divergence between the 
models of iLaw et al.l (120051 ) and iFellhauer et al.l (120061 ) in this region (see Figure [T9l middle 
panel), older material offers the strongest constraint on the model orbit. However, radial 
velocities for stars in the observed over-densities are required to confirm their association with 
Sgr. We are now obtaining radial velocities for RR Lyraes found along the debris stream. By 
confronting the radial velocity predictions of dynamical models with observations we hope 
to converge on increasingly accurate models for the Galactic halo and orbital elements for 
Sgr. 



8.3.2. Virgo Over Density 

The region between RA 180-230° (seen in Figure [201) exhibits two significant over densi- 
ties, one at RA=12'^47™, dec=-06°42" and distance of 16kpc and a second at RA=13''47™, 
dec=-ll°24™ and a distance of 19kpc. Both lie in the direction of several previously iden- 
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tifi ed halo substru c tures . The first indication of substructure in this reg ion was presented 
bv iNewberg et ali J2002I . region S297+63-20.5; at a distance of 18kpc). Izinn et al.l J2004h 
found a significant over density from their RRL survey (distance of 19kpc). Subsequent 
spectroscopy of a sample of RRLs and blue horizonta l branch stars frorn within this clump 
revealed a movi ng group wit h Vgsr = 100 ± 13kms~^ (IDuffau et al.ll2006l . dubbed the Virgo 
Stellar Stream). IJurig (120051 ) used photometric parallaxes of SDSS stars and found a broad 
over density centred on the above (the Virgo Over Density (VOD)) that covers ove r a thou- 
sand square degrees between distances of 5-15kpc. iMartinez-Delgado et al.l (120071) suggests 



that t he leading arm of tidal debris from the Sagittarius dwarf, as modelled by 



Law et al. 



( 2005 ) could provide the excess required to account for the VOD. Newberg fc Yanny ( 20051 ) 



have suggested the VOD may result from a non-axisymetric component to the Galaxy such 
as a triaxial halo. 



Newberg et al.l (120071 ) determines radial velocities for a sample of F stars drawn from 
the S297+63-20.5 region and finds radial velocities of Vgsr = 130±10kms~^. The low velocity 
dispersion of these stars is inconsistent with their origin in a triaxial halo and their radial 



velocity is at odds with that expected from stars in the Sgr leading tidal tail. iNewberg et al. 



( 120071 ) proposes that the three named features are one, we use the term Virgo Over Density 
hereafter to refer to this single system. 

The first oyer den sity, VOD Clump 1, lies 8 degrees southwest of the region identified 
by iDuffau et al.l (120061 ) as shown in Figure [5T1 This is followed by a second over density 16° 



further to t he south west, VOD C lump 2. It is apparent from the results of iBelokurov et al. 



(120061) and INewberg et al.l (120071 ) that regardless of the likely shape of the VOD the central 
regions must lie out of the SDSS region in the southern hemisphere. The twin clumps 
found by the present study plausibly form structure in the core of the VOD. However we 
must point out that due to the non-uniform survey coverage we are unable to clearly define 
spatial substructure. We are in the process of attaining radial velocities for samples from 



the tw o clumps to verify if they do indeed share the systemic velocity found by iDuffau et al. 
kood \ for the VOD. 



9. Summary 

In this paper we have presented a survey for RR Lyrae stars that has revealed a sample 
of 2016 RR Lyrae candidates from multiple (2-10) epoch data covering 1675 square degrees. 
Follow up photometric observations of these candidates show 76% to be RR Lyrae variables. 
Our survey reveals a series of significant over densities in the Galactic halo. Many of the 
over dense regions lie along the leading and trailing arms of debris from the Sagittarius 
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dwarf galaxy extending up to ~ 200° from the main body of the Sagittarius dwarf. Another 
locahty of over density can not be ascribed to the Sagittarius dwarf: namely the region of 
the Virgo Over Density. Here we find two significant over densities 10° to the south west of 
the centre previously reported. 

Radial velocities of RR Lyrae stars drawn from our over densities will give insight into 

the dynamics and origins of these systems. In particular, regions coincident with the leading 
and trailing arms of the Sagittarius dwarf could offer powerful constraints on the orbit and 
evolution of the Sagittarius dwarf. 

Wc thank Ken Freeman for his comments that helped to improve this paper. This data 
is based upon observations from Mount Stromlo Observatory's 1.3m Telescope before its un- 
timely demise in the bushfire of January 18th 2003. This paper utilises public domain data 
obtained by the MACHO Project, jointly funded by the US Department of Energy through 
the University of California, Lawrence Livermore National Laboratory under contract No. 
W-7405-Eng-48, by the National Science Foundation through the Center for Particle Astro- 
physics of the University of California under cooperative agreement AST-8809616, and by 
the Mount Stromlo and Siding Spring Observatory, part of the Austrahan National Univer- 
sity. This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 

Facilities: MtS:1.3m. 
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At (d) 

Fig. 1. — The time difference in days between epochs for the SEKBO fields. As expected 
from the survey design there are clusters of observations at 1-4 hours and 1-7 days apart. 
The group at 1 year are repeat observations of fields that failed the previous year. The insert 
shows a the distribution between 0.1- lOd. 
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Seeing (arcsec) 

Fig. 2. — The seeing distribution of images in the present study. The sohd hne is the 
cumulative distribution, while the dashed line shows the 3.6" seeing cut imposed on the 
observations as detailed in §4.21 
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Fig. 3. — Histogram of the number of images per field for the sample of 5212 fields. 




Fig. 4. — The distribution of fields in the present survey in galactic coordinates. The dashed 
line is the celestial equator. 
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Fig. 5. — The difference between the Landolt V (circles) a nd -R (crosses) standard magnitudes 
and those derived from our photometry for the IStetsoru (120001 ) stars in our sample. 
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Fig. 6. — The standard deviation of the photometry in Ryi as a function of magnitude for 
typical field in our survey. Overlaid is the sky-limited Poisson curve (line: see text) that is 
a fit to objects with unambiguously stellar shape parameters (the bold points). 
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Fig. 7. — The observed color variation of a sample of 119 MACHO RR Lyrae from the LMC. 
The line shows the best-fit relation AB = 1.32Ai? 
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Fig. 8. — RRL index L vers us V magnitude for one of the survey fields that overlaps the 
QUEST (IVivas fc Zinnll2006l ) survey. The two QUEST RRLs in the field are recovered above 
the dashed line that separates the non-RRLs of low L from the candidates with high L. 
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Fig. 9. — Typical completeness profiles as a function of V magnitude for fields with 2, 3, 4, 
and 5 epochs (a to d respectively). The shaded histograms represent RRL type cs. 
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Fig. 10. — A map of completeness across the survey. The color represents the V magnitude 
at which the completeness drops to 25% in the case of the RRab (top two panels) and RRc 
(bottom two panels) variables. The celestial equator is shown as the inclined dashed hne. 
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Fig. 11. — The asymmetric nature of the RRLa6 hght curve introduces a magnitude depen- 
dent offset to the mean magnitude derived for our RRL candidates. The mean magnitude 
offsets for Bm (the sohd hnc for fields with two observations; the dotted hne for more than 
two observations) and Ru (the long dashed line for two observations; short dashed for more 
than two) are shown as a function of magnitude for the ensemble of RRL candidates. 
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Fig. 12. — The classification of 66 RRL candidates obtained from follow-up photometry. 
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Fig. 13. — (a) The normalised distribution of candidate RRL colors compared to that of 
a representative sample of ~11000 dereddened b > 60° halo stars, (b) the corresponding 
color-magnitude diagram for the candidate sample and halo population. The long dashed 
hne is the mean V—R of the MACHO RRL distribution discussed in §5. The dotted hne at 
V—R—0.3 is the imposed color cut. 
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Fig. 14. — Left: The radial distribution of RRL candidates in ecliptic longitude. Right: The 
corresponding distribution of objects from one of our monte-carlo simulations of the halo (see 
§8.3p . A number of substructures (shaded in the figure) are apparent in the observations that 
are not found in the simulation, in particular at ecliptic latitude 120°, 190°, 210°, and 320°. 
The significance of these structures is discussed in §8.3[ The circles are placed at ^=15, 17 
and 19. 
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Fig. 15. — Our 2016 RRL candidates color coded by distance overlaid upon the view of 
the outer halo due to iBelokurov et al.l (120071 ) which includes analysis of th e main-sequence 



turnoff stars from SDSS (the top right insert) and M giants from the study of lMajewski et al. 



(120031 ) (the greyscale background). The Sagittatius Stream dominates the view of the outer 
halo (the main body of the Sagittarius dwarf is located centrally in this figure at Dec~ —20°. 
See the electronic edition of the Journal for a color version of this figure. 
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Fig. 16. — Variation of the fractional halo volume as a function of Galactocentric distance, 
f{R). Left panel: for two fields centred on {£,b) = (213°, 34°) (solid line) and (262°, 60°) 
(dotted line). The asymptotic limit represents the constant solid angle subtended by the 
fields if i?o = with infinite limiting magnitude. Right panel: The same profiles corrected 
for completeness as a function of R. The profiles for all fields are summed to give the total 
f{R) for the survey. 
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Fig. 17.- 



Total fractional volume of the halo surveyed by the 3692 fields 
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Fig. 18. — The apparent spatial density of RRL candidates determined assuming that all 
candidates are RRLs (see text for discussion of the best- fit lines). We use this result to 
quantify the average background upon which substructure resides. 
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Fig. 19. — Top: A significance map of over densities in RRLs in tlie lialo (see text for 
details on derivation). Regions in black represent over densities of 80" significance. The 
circles are at heliocentric distances of 20, 40, 60 and 80kpc. The dashed straight line is 
the line of nodes for the intersection between the orbital plane of the Sagittarius dwarf 
galaxy (llbata et al.ll2003l ) and the echptic along which our survey is performed. The dotted 
lines are distances at which the plane of the Sagittarius dwarf is 1.5, 3, and 6kpc away 
from the survey volume. Mid dle: A view of the Sagi ttarius plane showing M giants from 



Majewski et al.l (|2003l. points) . iBelokurov et al.l (l2006l . squares) distances and de tections of 



Sgr st ars from iNewberg et al 
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Fig. 20. — Close-up view of the inner 40kpc of our survey. While much of the structure can 
be attributed to the debris of Sgr, the Virgo Over Density region at RA~ 180°-210° and 
distance ~ 20kpc is evident. Circles are at heliocentric distances of 20 and 35kpc and angles 
are in degrees of right ascension. 
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Fig. 21. — The distribution of RRL candidates in the region of the Virgo Over Density 
covered by our survey (survey fields are the open squares). The points shown have heliocen- 
tric distances of 16 < r < 26kpc. Clump 1 is seen at RA=192°, dec=-8° an d Clump 2 at 



RA= 206°, dec=-12°. The sample of RRLs and blue horizontal branch stars of iDuffau et al. 



( 120061 ) are shown as filled circles and triangles respectively. The dashed outline shows the 
region where Duffau et al. found an excess of main-sequence stars of appropriate apparent 
magnitude to be associated with the Virgo Over Density. 
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Table 1: Wide-field Surveys for RRLs. 



Survey 


Sky Coverage 


Number of 


Maximum 


Name 


(sq. deg.) 


Candidates 


Distance(kpc) 


QUEST 


380 


498 


55 


LONEOS-1 


1430 


838 


30 


Ivezic et al. (2000) 


100 


296 


65 


Sesar et al. (2007) 


290 


634 


100 


Current Survey 


1675 


2016 


50 



-48- 



Table 2: Table of RR Lyrae candidates found in the present work. 



Object ID 


RA 


Dec 


V 


V-R 


E(V-R) 


Heliocentric Distance (kpc) 


128416.544 


00:12:52.83 


05:02:56.71 


17.17 


0.19 


0.014 


20.3 


116375.115 


00:13:28.83 


02:03:43.68 


14.50 


0.08 


0.015 


5.9 


128412.1045 


00:14:37.09 


03:51:27.99 


18.35 


0.20 


0.013 


35.0 


115608.149 


00:16:14.73 


01:53:51.58 


15.34 


0.18 


0.019 


8.6 



Note. — Table [2] is published in its entirety in the electronic edition of the Astrophysical Journal. A 
portion is shown here for guidance regarding its form and content. 



